Multiscale approach to spin transport in magnetic multilayers. 
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This article discusses two dual approaches to spin transport in magnetic multilayers: a direct, 
purely quantum, approach based on a Tight-Binding model (TB) and a semiclassical approach 
(Continuous Random Matrix Theory, CRMT). The combination of both approaches provides a 
systematic way to perform multi-scales simulations of systems that contain relevant physics at scales 
larger (spin accumulation, spin diffusion...) and smaller (specular reflexions, tunneling...) than the 
elastic mean free paths of the layers. We show explicitly that CRMT and TB give consistent results 
in their common domain of applicability. 

PACS numbers: 72.25.Ba, 75.47.-m, 75.70.Cn, 85.75.-d 



I. INTRODUCTION 

Many interesting physical effects can be observed in 
non collinear magnetic systems where the various mag- 
netizations are not aligned or even not in the same plane. 
Indeed, in addition to the Giant Magneto Resistance 
(GMR 1 - 2 ) and Tunneling Magneto Resistance (TMR 3-5 ) 
effects observed in the collinear configuration, the non 
conservation of spin current gives rise to an influence 
of electronic transport on the magnetization dynamics 
which itself can lead to various phenomena such as mag- 
netization reversal 6-9 or dc driven radio frequency os- 
cillators 10-14 . The potential for applications opened by 
these effects (magnetic memories, reprogrammable logic, 
tunable rf sources...) has been recognized very early and 
lead to an important expansion of the field 15 . 

Consequently, the theory of non collinear magnetic sys- 
tems has received a lot of attention. Many different ap- 
proaches have been developed, ranging from purely quan- 



tum ' to semi-classical Boltzmann equation 
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(gen- 



eralized) circuit theory 2 1-27 or Random Matrix Theory 
(RMT) 28 ' 29 or full ab-initio (density functional theory) 
calculations 30-33 . A good understanding has now been 
reached of the respective domains of applicability and 
links between the various approaches 26,29 . On the other 
hand it is becoming increasingly clear that the magneti- 
zation dynamics of real devices cannot be properly cap- 
tured by simple (analytically tractable) models and that 
numerical simulations that treat the magnetic and trans- 
port degrees of freedom on an equal footing need to be 
developed. Several steps have already been taken in that 
direction 34-38 and should eventually lead to simulations 
of good predictive capabilities. Let us mention, as ex- 
ample, the recent experiments on spin torque induced 
ferromagnetic resonance in the nonlinear regime 39 which 
has been successfully compared to coupled simulations 
at the macrospin level for the magnetic part and one- 
dimensional semi-classical level for spin transport. Other 
geometries with stronger spin texture (domain wall, vor- 
tex, strong oersted field...) will require a full three di- 
mensional treatment however. 



This paper is devoted to a discussion of a multiscale 
approach to the spin transport theory of magnetic mul- 
tilayers: for some system, such as magnetic tunnel junc- 
tions, semi-classical approaches fail (they are intrinsically 
Ohmic and cannot capture the exponentially small trans- 
parency due to the presence of the insulating oxides) and 
one is therefore tempted to use fully quantum mechanical 
simulations. Those simulations are extremely computa- 
tionally costly however, and cannot be performed on full 
size realistic devices. The multi-scale approach taken in 
this paper consists of describing the intrinsically quantum 
parts of the system at the quantum level (with the help of 
a Tight-Binding TB model) and embedding those parts 
in a semi-classical description (valid on scales larger than 
the elastic mean free path). This multi-scale approach is 
performed in three steps. 

First, we provide an extensive presentation of our semi- 
classical approach which we refer as Continuous Random 
Matrix Theory (CRMT). CRMT was developed recently 
in Ref 29 to which we refer for a compact account of the 
theory. We discuss in details the links with other ap- 
proaches (in particular the explicit mappings to circuit 
theory 26 and Valet- Fert theory 40 ) and explain how clas- 
sical concepts like spin accumulation naturally appear in 
the formalism, even in the absence of local equilibrium. 

Second, we study a toy model for a tunneling spin valve 
and show explicitly how to combine CRMT with the TB 
approach. In particular, we discuss in details the condi- 
tions for the CRMT approach to be valid by comparing 
the multi-scale approach with full quantum calculations 
of the TB model. Third and last we construct an effec- 
tive TB model capable of describing diffusive magnetic 
multilayers. This TB model is solved at a purely quan- 
tum level and we find a close agreement between TB and 
CRMT, so that the same parametrization can be used 
for both. The appendix describes a direct integration 
of CRMT equations for a magnetic layer with no spin 
texture, together with a discussion of the boundary con- 
ditions and role of Sharvin resistances. 
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II. CONTINUOUS RANDOM MATRIX 
THEORY (CRMT) 

This section is devoted to the CRMT semi-classical 
theory that describes the conducting electrons of a mag- 
netic multilayer. This section can be considered as a 
long version of the letter published in Ref 29 to which 
we add an intuitive derivation of the relation between 
Scattering matrices and spin accumulation, making use 
of the voltage probe concept 41-43 . CRMT is a continu- 
ous version of a Random Matrix Theory (RMT 44 ) for 
non-collinear magnetic multilayers that was developed 
in Ref 28, which is itself based on the fully quantum 
Landaucr-Buttiker scattering approach 41 . For the sake 
of completeness, both are briefly reviewed below. 

A. Landauer Buttiker approach to magnetic 
multilayers 

The Landauer-Buttiker or scattering matrix 
formalism 45-47 is a standard approach to quantum 
transport. A sample is defined by the scattering matrix 
S which expresses the outgoing propagating modes in 
term of the incoming ones. The incoming states are 
filled according to the Fermi-Dirac distribution of the 
electrodes to which they are connected which enables to 
relate various physical quantities, such as conductance 
G or spin current J, to the matrix elements of S. 
The system contains N c h 3> 1 propagating modes per 
spin where _/V c h ~ A/\ 2 F (A is the transverse area of 
the electrode and Xp is the Fermi wavelength). The 
amplitude of the wave-function on the different modes is 
given by a vector ipi±, 

where ipi±a is a iV c h vector that contains the amplitudes 
for right (left) moving electron direction with spin a =t, -J, 
along the z-axis in region i = 0, 2, see Fig. 1 for a cartoon. 
With these definitions, the S matrix is defined as, 

The 5* matrix is a 4-/V c h x 4-/V c h unitary matrix and con- 
sists of 2A^ c h x 27V c h transmission t, t' and reflection r, r' 
sub-blocks, 




while the transmission and reflection matrices have an 
internal spin structure: 

-Gsfc) (4 » 

where t aa > are N<± x N<± matrices containing amplitudes 
for transmission between u' and a spin states. 




FIG. 1: Cartoon of the S matrix approach. We define the 
region on the left and 2 on the right (region 1 in between 
will appear later in the text) . The + (-) respectively stands for 
right (left) going modes. The mode amplitudes are denoted 
tpi± a at the quantum level and become probabilities Pi± a in 
RMT. 

The conductance of the system is given by the standard 
Landauer formula G = ^-Tr [th] while the spin current 
in region is given by, 

where a is the vector of Pauli matrices a = {S x ,a y ,a z ) 
and [i the difference of chemical potential between the 
two electrodes (with similar formulas for the spin cur- 
rent inside the system 28 ). Note that in addition to the 
non-equilibrium spin current, there can exist an equilib- 
rium one which in the context of magnetic multilayers is 
often referred as the interlayer exchange (RKKY type) 
coupling 48-51 . 

B. Random Matrix Theory (RMT) for magnetic 
multilayers 

So far, the theory is fully quantum and contains in 
particular many interference effects such as weak localiza- 
tion 52 ' 53 or universal conductance fluctuations 54 ' 55 . We 
now proceed with the description of the semi-classical 
RMT 28 theory. We give a presentation slightly more gen- 
eral and compact than in Ref. 28 to which we refer for 
the proofs. The systems in which we are interested in 
this paper, typically pillars of a few tens nanometers in 
diameter connected to top and bottom electrodes, con- 
tain many channels, typically Af c h w 10 4 — 10 5 . When 
the scattering is not perfectly ballistic (Fermi momentum 
mismatch at the interfaces, surface roughness or impurity 
scattering) those channels get mixed up. RMT assumes 
that this mixing is ergodic 44 , (i.e an electron entering the 
system in a given mode will leave it in an arbitrary mode 
and pickup a random phase in the process). The the- 
ory can be derived systematically when the system size 
is larger than the mean free path, but we shall see that 
in practice, its domain of applicability is even wider. 

RMT is expressed in term of 4 x 4 " hat" matrices that 
are obtained by tracing out the transverse degrees of free- 
dom of the original reflection and transmission matrices. 
For instance, t is defined as, 

i<rq,<r'7\' = ~T T Trjy ch \taa'tL ,\, (6) 
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(TrAr ch trace is taken on the modes only) or equivalently, 
giving the 4x4 structure explicitly, 



t = 



N. 



-Tr 



ch 



( ttt4t *tt*tj. *t-l-*|t *t4-*|j. \ 

*tt*it *tt*u *n*Lt hi^u 

http htt+i *u*n t n t n 

V ^t f it *+T*U *-W-*4-t / 



(7) 



Similarly, "hat" -matrix S 1 has a form similar to Eq. (3), 

In analogy with "hat" matrices, we introduce 4- vectors 
Pi± which corresponds to the modes amplitudes r 4>i± IJ : 



Pi± = 




(9) 



The components of the 4-vector Pi± have interpretation 
in term of probabilities. For instance, Po+t (P2-4.) ac- 
counts for the probability to find a right (left) moving 
electron in region (2) with spin \ (1). The mixing com- 
ponents, P mx are complex numbers which correspond to 
probability to find the electron along the x (real part) or 
y (imaginary part) axis. Inside magnetic layers where the 
z axis will correspond to the direction of the magnetiza- 
tion, they will correspond to the (small) probability for 
the spin to have a part transverse to the magnetization. 
Again, in analogy with Eq. (3) which expresses the am- 
plitudes of the outgoing modes in term of the incoming 
ones, we have, 



Po- 
P+ 



S 



Ph 
P 2 



(10) 



When the " mixing" components of the 4- vectors play no 
role, as in collinear systems, Eq. 10 has an obvious in- 
terpretation in term of a Master equation. For instance, 
its first row expresses that the probability to find a left 
going electron in region has two contributions coming 
from the probability to have a reflection and transmission 
event. In other word, instead of the original interference 
problem with amplitudes, one now deals with the classi- 
cal equivalent with probabilities. For non-collinear sys- 
tem however, the presence of the quantum SU(2) struc- 
ture of the spin introduces some complex numbers (the 
mixing coefficients) and the Markov process interpreta- 
tion does not hold, strictly speaking. 

Physical observables. The chief result of Ref. 28 are the 
expression for the currents and spin current in term of the 
"hat" matrices, to leading order in N c h. Equivalently, we 
can write those expression with the P-vectors, 



Ji 



47T 



a P 



i+ 



(11) 



and 



I = 



+t 



+1 



P-l] 



(12) 



where I is the charge current and lZ s h = h/(N c i 1 e 2 ) 
is the Sharvin resistance. <y = (<ff\-, a^, a^i, an) is 4- 
vector composed of components of Pauli matrices. To 
complete the theory, we need the boundary conditions 
imposed on the incoming electrons on both sides of the 
system. We focus here on normal electrodes, but these 
conditions can be easily generalized to the case where the 
electrodes are themselves magnetic. We have, 



0+ 





(13) 



where /Uo and ^2 are the respective chemical potentials of 
the two electrodes. For a given "hat"-matrix S, the com- 
binations of Eqs. (10,11,12) and Eq. (13) form a complete 
set of equations to obtain the physical quantities. 

Two systems in series. Usually, the S matrix can be 
obtained for parts of the system (interfaces or the bulk 
part of one material) and it is important to be able to 
combine several parts to obtain the global S matrix. This 
is particularly important when one is interested in calcu- 
lating the spin torque, as one also needs to calculate the 
spin current flowing inside the system and not only near 
the electrodes. Let us consider the transport through 
a system which consists of two subparts, as pictured in 
Fig. 2. The transport through each subpart is charac- 
terized by the corresponding "hat" -matrix S a and Sb- 
Eq. (10) reads for each of them, 



Po- 
Pi+ 



S a 



0+ 



Pi 



Pi- J ' V p 2+ ) - Sb V P2 



1+ 



(14) 



Eliminating Pi we find addition law for "hat" -matrices: 



Po- 

P2+ 



— S a +b 1 p 



0+ 



1 - r b r a 



fa+b = f b + t' b 



1 - f a r' b 



with, 



r a tb 



(15) 
(16) 
(17) 



and similar expressions for f' a+b and i' a+b . We can also 
use Eq. (14) to express Pi± in terms of the incoming 
fluxes: 



Pn 
Pi 



1 P I f P 

l-r a f' b b a l-f a f' b ' aib 
1 ftp 1 P 

i-f'„fa b a i-f' b fa b 




(18) 



Using Eq. (11), Eq. (13) and Eq. (18), one can calculate 
the spin current in the region 1 inside the system. 
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FIG. 2: 

Cartoon of a system made of two subsystems a and b. 

Rotation matrices. The scattering matrix of a piece 
of magnetic material is most easily found in the ba- 
sis parallel to the local magnetization, i.e. the reflec- 
tion and transmission matrices are given for the major- 
ity and minority electrons. For non-collinear multilay- 
ers where different magnetization directions come into 
play, one needs to rotate the original S matrix onto its 
form S = R$fiSRg - in the chosen working basis. Here 
Re.n = cxp(— iff-n 9/2) — cos(9/2) — ia-n sin(0/2) is the 
rotation matrix of angle 9 around the unit vector n that 
brings the magnetization onto the z-axis of the working 
basis. In term of "hat" matrices, this translate directly 

into S = RgfiSRl - with, 

Rarj.a'rj' = R<j<j' R^rf i (1^) 

a unitary matrix. 

C. From scattering degrees of freedom to spin 
accumulation and spin current 

The natural variables of RMT as it was introduced 
above are the 4-vectors P x ± which characterize the 
"probability" to find a left or right moving electron in 
region x. Let us now introduce a new set of variables 
defined as, 

j(x) = [P + (x)-P-(x)}/(e1l sh ), (20) 

/z(z) = [P+(z)+P-(z)]/2 (21) 

As we shall see, these two 4-vectors correspond respec- 
tively to the (spin resolved) current and chemical poten- 
tials flowing in the system. There are two complementary 
ways to make this connection. The first one is to write 
our fundamental equations (essentially Eq. (10)) in term 
of these new variables. This will be done toward the end 
of this section, and we will find that j(x) and fi(x) sat- 
isfies the well known Valet-Fert diffusive equations 40 for 
continuous collinear systems and the equations of circuit 
theory for discrete non-collinear ones 24 . This connec- 
tion is very interesting from the theoretical point of view 
as very different routes have been taken to obtain these 
equations (scattering and Random matrix theory on one 
side and Keldysh Green function formalism and quasi- 
classical approximation on the other side). It has also a 
practical interest as one can use Eq. (20) and Eq. (21) to 
go back and forth between a " scattering" approach and a 



" diffusive" approach and both have technical advantages 
for practical calculations. 

There is a second, more direct, way to connect respec- 
tively j(x) and n(x) to the concepts of spin currents and 
spin accumulation. For j(x), the connection is straight- 
forward, as Eq. (11) reads, 

f{x) = lff-j(x) (22) 

while 

i=h+h ( 23 ) 

so that j-c and jj. can be directly interpreted as spin cur- 
rents for up and down electrons (while j mx accounts or 
spin current transverse to the z-axis). The connection 
between n(x) and a hypothetical spin resolved chemical 
potential is more problematic as scattering theory does 
not have any notion of local chemical potential. In fact, 
the existence of a local chemical potential would imply 
some sort of local equilibrium inside the system. How- 
ever, the theory that we have developed so far is purely 
elastic and such a local equilibrium is not present. We 
will find that everything happens as if there were some 
sort of local equilibrium, except for one important point: 
the presence or absence of the contact resistance. 

In order to provide a physical understanding of ji{x), 
we make use of a theoretical trick known as the volt- 
age probe 41-43 : we connect a point inside the system 
(region 1) to an external electrode and adjust the (spin 
resolved) chemical potential of this electrode so that no 
(spin) current flows from/to it. In the absence of mag- 
netism, this could be achieved experimentally with the 
help of a STM tip. In mesoscopic physics, the voltage 
probe is used to induce some decoherence 56 in an other- 
wise perfectly coherent theory, hence introducing a finite 
phase coherence time which depends on the coupling of 
the system with the probe. Here, we will take this cou- 
pling to be extremely small, so that the probe will not 
affect the physics of the system. On the other hand, we 
are interested in the voltage that one must apply on the 
probe to stop any current from flowing from/to it: this 
will be our definition of an effective local chemical po- 
tential, and we shall find that this definition matches the 
one of n{x). The setup is shown in Fig. 3. In addition 
to S a and Sb defined previously, we introduce the matrix 
S e which describes the (small) coupling to the probe, 

/ (l-e)l el \ 
S £ = (l-e)i el , (24) 

V ei el (1 - 2e)i / 

with eCl. When e = the probe is decoupled from the 
conductor and the electrons freely propagate between the 
regions 1 and V . The full system is entirely determined 
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i+ 




2+ 



FIG. 3: (a) Cartoon of the system in presence of the voltage 
probe. 



by the following set of equations, 



P _ \ _ A ( Ph 
P 1+ - ba { P 




P 2+ J - Sb ( iV 



(25) 



to which one must add the incoming boundary conditions 
on the three electrodes Po+, P2- [Eq. (13)] and P p + (the 
probe). Eliminating P\ t y we obtain S. In particular we 
get to leading order in e, 



S p ,o = e(l + f' b ) 
Sp,2 = e(i + f a ) 



1 



1 - f a f' b 
1 

i - f' b f a 



1 -2el + 0(e 2 ) 



(26) 



We now impose that J p = [Eq. (11)] which provides 
our boundary condition on the probe: 



Using Eq. (18) and Eq. (21), we arrive at, 
Pp + = \(Pi++Pi-) = »i 



(27) 



(28) 



In other words, the (spin resolved) voltages that one need 
to apply on the probe is equal (as announced above) to 
the formal chemical potential that we have defined in 
Eq. (21) so that it is legitimate to call fi(x) a chemical 
potential. 



interfaces between several magnetic and non-magnetic 
layers 57 . In the same spirit, in the next section, we will 
calculate the "hat" matrices from calculations of an effec- 
tive quantum model for diffusive magnetic metals. In the 
second route, that we take here, we derive those "hat" 
matrices from phenomenological considerations. The S 
matrix for a piece of bulk material or for an interface 
between two different metals will be parametrized by a 
few parameters. These parameters will in turn be put 
in one to one correspondence with the parameters of the 
(well established) Valet-Fert theory 40 . This has a double 
advantage as it allows to make a direct connection with 
a widely spread theory, and it also allows to make direct 
use of the important experimental effort that has been 
done to parametrize Valet-Fert theory. 

General form of the "hat" matrices in a magnetic 
material. A typical magnetic multilayer consists of al- 
ternating layers of magnetic and non-magnetic materi- 
als of various widths. We decompose the correspond- 
ing S matrices in the bulk and interface parts that will 
be parametrized independently. The general form of a 
"hat" matrix (t in what follows but similar considera- 
tions apply to t' , r and r') is a full 4x4 matrix given by 
Eq. (7). In the absence of spin-orbit scattering (and/or 
magnetic impurities), spin is a good quantum number of 
the problem and the t matrix is diagonal in spin space 
(in the basis parallel with the magnetization). Hence t 
is also diagonal. It consists on one hand of the prob- 
abilities Tff (Tn) f° r an t (■!•) s P m t° be transmitted 
and on the other hand of the so-called "mixing trans- 
mission", T mx = (l/7V ch )Tr Nch (< n t^). The latter is a 
complex number. Its amplitude measures how much of 
a spin transverse to the magnetic layer can be trans- 
mitted through the system while its phase amounts for 
the corresponding precession 58 . Due to large number of 
modes the average T mx decays rapidly with the size of 
a ferromagnet. 5 ' 58 Thus the mixing terms are small in 
magnetic systems, but can (and do) play a role in non- 
collinear configurations nevertheless. Spin-orbit interac- 
tion however cannot be entirely neglected, as it leads to a 
finite spin diffusion length in the system. Hence, a finite 
probability T^j. (T^) for a 4- (t) spin to be transmitted 
as a t (4) spin must be considered. In the following, we 
suppose that mixing elements (complex numbers) that 
also involve spin-flip scattering can be totally ignored, 
leading to the following form of the "hat" -matrix: 



D. Continuous Random Matrix Theory (CRMT) 

The basic ingredient of the RMT theory presented 
above are the "hat" scattering matrices of the various 
sub parts of the system. Once those are known, they 
can be concatenated to obtain various physical observ- 
ables. There are two routes that one can take to obtain 
those "hat" matrix. The first one is to go back to their 
definition in term of the original quantum problem, and 
compute them from purely quantum calculations. This 
route has been taken with ab-initio calculation for the 



/ T n 








T n 



















mx 





w 








Tu 



(29) 



with the probabilities T aa ' = l/i\ r c hTrjv ch (^. a ./t ( r<r')- Note 
that the unitarity of the S matrix imposes the following 
constraint on the mixing coefficient, 



\T mx \ < \Att^44 



(30) 



General form of the "hat" matrices in a non-magnetic 
material. In a non-magnetic material there is no pre- 
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ferred direction for the spin thus a "hat" -matrix of the 
non-magnetic metal must be invariant with respect to ro- 
tations, as defined in Eq. (19). One finds by inspection 
that there are two 4x4 matrices invariant with respect 
to rotation around an arbitrary axis: 



r t , r 4 r sf and r E 



h = 



( 1 

10 

10 

\0 1 



10 1 





10 1 



(31) 



Hence, we write the "hat" matrix of a normal layer as a 
combination of these two invariants: 



i=[T-2T si ]h+T s{ i 2 . 



(32) 



where T is the total transmission probability while T s { is 
the probability of transmission with spin-flip. 

Obtaining S for bulk materials. So far we have dis- 
cussed the form of the "hat" matrices for an arbitrary 
subpart of the system. This general form is directly use- 
ful for the parametrization of the "hat" matrices of in- 
terfaces. For bulk material, we introduce the S matrix 
of a very thin slice of material of width SL. S(SL) (and 
therefore the bulk properties of the material) is entirely 
characterized by two matrices A* and A r defined as, 



t(6L) = 1 - A l 5L , f(6L) = A r 6L 



(33) 



Once S(SL) is known, one can make use of the addition 
law derived in Eq. (16) and Eq. (17) to obtain a dif- 
ferential equation that allows the computation of S(L): 
combining S(L) and S(5L) one obtains S(L + SL) which 
by taking the limit SL — > provides, 



df 
dL 
dt_ 
dL 



= A r - A*f - rA* + fA r f (34) 
= -AH + rA r t (35) 



Eq. (34) and Eq. (35) can be integrated to obtain S(L) 
of an arbitrary bulk part (see the Appendix). However, 
for numerical purposes, it is more efficient to use directly 
Eq. (16) and Eq. (17) which leads to an extremely fast 
integration time cx log L (In practice, one starts with an 
extremely small piece of material described by Eq.(33) 
and recursively doubles it size using Eqs. (16,17) until 
obtaining the full width L of the layer) . 

The matrices A* and A r are parametrized by the gen- 
eral form given in Eq. (29) an Eq. (32). The parametriza- 
tion is further constrained by the unitarity of S (conser- 
vation of current). Due to the origin essentially ballistic 
of the finite mixing coefficient, we neglect the correspond- 
ing contribution to A r . Eventually, a bulk magnetic ma- 
terial is characterized by four independent parameters 



A* = 



A r = 



/r t + r sf o o -r sf 

o r mx o o 
o o o 

V -r sf o o r+ + r sf 

/r t -r sf o o r sf \ 





V r sf oor r r sf J 



(36) 



(37) 



These four parameters correspond in turn to 5 different 
lengths. The two most important one are the mean free 
paths for majority (If) and minority (Z4.) electrons de- 
fined as l a = l/r<j. Next comes the spin diffusion length 
l si = [Ar^ir^ + r^]- 1 / 2 (see below Eq. (46) and the Ap- 
pendix). Last the complex number T mx = l/l± + i/lh 
where Z j_ is the penetration length of transverse spin cur- 
rent inside the magnet while Zl is the Larmor precession 
length. These definitions are supported by the form of 
T mx (L) which is readily obtained by integrating Eq. (35), 



T mx (L) = e 



_ p -L/l ± -iL/l^ 



(38) 



l± and Zl, are believed to be roughly equal, and the small- 
est characteristic lengths with typical values in the nm 
range. The explicit form of T acr > (L) will be discussed in 
the appendix. 

In a normal metal, the parametrization obeys Eq. (32) 
and we have 



a* = [r + 2r sf ]/ 1 -r sf 7 2 . 



A r = [r-2r sf ]/ 1 + r sf / 2 . 



(39) 



(40) 



The theory is now formally complete. A given multi- 
layer is then constructed by using the addition law Eqs. 
(16,17) for the various bulk layers and the corresponding 
interfaces. 



E. Correspondence between CRMT and Valet Fert 
theory 

We are now ready to connect CRMT with other ap- 
proaches taken in the literature. An important connec- 
tion that we make here is to show that for collinear sys- 
tems, CRMT equations are equivalent to the celebrated 
Valet-Fert (VF) equations. 40 As these equations have 
been extensively parametrized with the huge corpus of 
experimental data available in CPP GMR, this allows us 
to transfer directly this parametrization to CRMT. 

Bulk magnetic materials. Let us start with S(Sx) 
which relates P±(x) and P±(x + Sx) on the two sides 
of a thin slice of material according to Eq. (10). The 
explicit form of S(6x) given by Eq. (33) provides differ- 
ential equation for P±(x) which are the counterparts of 
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Eq. (34) and Eq. (35): 
dP+(x) 



dx 
dP-(x) 
dx 



= -A t P+(x)+A r P^(x) 
= -A r P+(x)+A t P^(x). (41) 



In the absence of magnetic texture, the equations for P±,t 
and P±.i are decoupled from the equations for P± im x and 
P± mx . Focussing on the non-mixing contributions (first 
and last row of the 4- vectors) of Eq. (41) we get explicitly, 



d_ 

dx 



= -A* 



P- 



i i (42) 



where A* and A r are 2x2 matrices, 



A t/r = 



r t ± r 8f 
Tr sf 



Ti ± r sf 



(43) 



Eq. (42) accounts for the conservation of probability in 
the scattering events, i.e. it is the Master equation of the 
underlying Brownian motion undertaken by the electrons 
upon the various reflection and transmission events. 

Using Eq. (20) and Eq. (21) we can now write Eq. (42) 
in terms of j(x) and n(x), and arrive at, 

= -l/(eT a 1Z sh )d xf i m (x) (44) 
9x3t/d x ) = 4r sf /(eft sh ) K/tW - f^/i( x )] ( 45 ) 

which are precisely the Valet-Fert equations 40 . Hence, 
for a collinear system, CRMT simply reduces to VF the- 
ory. This allows us to build a one to one correspondence 
between the CRMT parameters (Tf, T± and r s f) and 
the VF parameters \pf, pi (resistivities for majority and 
minority electrons) and l s f (spin-flip diffusion length)]. 
Using the standard notations for the average resistivity 
p* and polarization (3 \p^(\) — 2p*(l =F P)] we find 



1 



T^sh 




(46) 
(47) 
(48) 



Note that the mixing coefficient r mx is not fixed by this 
parametrization as they only play a role in non collinear 
configurations. So far there are very few experimental 
data 59-61 allowing to extract r mx so that one often relies 
on model or ab-initio calculations to estimate it. 

Interfaces. In VF theory, the interfaces are charac- 
terized by effective interface resistances r-j- (rj,) for ma- 
jority (minority) electrons as well as the spin-flip prob- 
ability S. Alternatively one can introduce the average 
resistance r b * and polarization 7 [r^j, = 2r b *(l =p 7)]. 
Much is known experimentally 62,63 about r b j_, but there 

are much less experimental data 64 for S (and those are 
mainly for normal- normal interfaces). In VF theory, the 



interface boundary conditions arc obtained by introduc- 
ing an virtual material of width d, resistivity p e J' and 
spin diffusion length ■ This virtual material is then 
taken to be infinitely thin d — > while keeping the in- 
terface parameters finite: 6 = d/l e s (* and r b a — p e J? d. 
Repeating the same procedure for CRMT allows us to 
map the VF parameters to CRMT interfaces, 



(l + e- s )/2 



1 + 2(r b * /K sh ){l - 7) 

(l-^)/2 
l + 2(r & */^)(i~7) 

(1 - e~ s )/2 
l + 2(r b */K sh )(l + 1 ) 

(l + Q/2 

l + 2(r 6 */^)(l + 7) 
1 



= 1 - 



(49) 
(50) 
(51) 
(52) 

,^t = (53) 



\ + 2{r h */n sh ){i- 1 y 
u \ . = 1 - i 1 Knini -1 , ~ ,Ru = (54) 



\ + 2{r b */n sh ){i + 1 y 



(55) 



Alternatively, the interface S matrix can be obtained di- 
rectly from ab-initio or model quantum calculations (see 
the next section). Once again, the mixing reflection and 
transmission coefficients are not fixed by CPP GMR ex- 
periments. Those numbers can be extracted from angu- 
lar resolved magneto-resistance or spin pumping experi- 
ments but limited data are available so far 25 . 

Valet-Fert theory and Sharvin resistances. We have 
just proved that for a collinear system, VF theory and 
CRMT are simply equivalent. There is however a small 
difference that appears in the boundary conditions at the 
electrodes. In CRMT those boundary conditions come 
from the (quantum) Landauer formula and the presence 
of different voltages between the reservoirs located at, 
say, x = and x = L imposes P +a (0) = eV(0) and 
P- a (L) = eV(L) [sec Eq. (13)]. A direct consequence 
of these boundary conditions is the existence of a fi- 
nite resistance, even for perfectly transparent interfaces 
and materials with negligible resistivity. In the context 
of mesoscopic physics, this leads to the quantization of 
conductance in unit of 2e 2 /h which has been observed 
repeatedly. 65 ' 66 . These boundary conditions can be ex- 
pressed in term of p a (x) and ja(x), and give, 



MO) 



■(67^/2)^(0) 
(eK sh /2) j a {L) 



eV(0) 
eV(L) 



(56) 
(57) 



In other words, one needs to add 1Z s h/2 resistors on the 
two sides of the multilayer. For typical spin valve pillars, 
the intrinsic resistance of the pillar is only a few time lZ a h 
so that one really needs to take into account the presence 
of these Sharvin resistances in series to properly describe 
the pillar. We note that those proper boundary condi- 
tions can easily be included in a standard VF calculation. 
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A detailed discussion of the effect of 7?. s h on the polariza- 
tion of spin current will be discussed in the Appendix. 

F. Correspondence between RMT and 
(generalized) Circuit Theory 

Let us now turn to non-collinear configurations. An 
alternative popular and powerful approach in this case 
is the so-called circuit theory 24 where the system is dis- 
cretized into various parts connected by "nodes" where 
one defines the spin resolved chemical potential. Circuit 
theory, initially derived for very resistive elements has 
been further extended into the "generalized" circuit the- 
ory to properly take into account the Sharvin resistance. 
The similitude between (generalized) circuit theory and 
RMT was recognized very early and it was shown in many 
cases that both theories gave the same result 26 . Here we 
show that the analogy is in fact complete: RMT and 
generalized circuit theory are the same theory written in 
different variables: P + and P for RMT and /x and j 
for circuit theory. One simply goes from one to the other 
using Eq. (20) and Eq. (21). 

The most general version of RMT between the two 
sides of a conductor is given by Eq. (9). Turning now to 
/x and j variables on the two sides L (left) and R (right) 
of the conductor, we get, 

(1 + f + oft)-^-3 R = 26' n L - (1 - f + 6't)n R (58) 

(1 + f' - ot')^h = (l-f' + 6t')» L - 2o» R (59) 

where 6 = t(l + r) -1 and & = t'(l + f') _1 - These two 
equations can be seen as the generalization of Ohm law 
and their linear combination provides the conservation 
equation for spin current. Eq. (58) and Eq. (59) can be 
considered as an extension of (generalized) circuit theory, 
including in particular interface spin-flip scattering. 

To recover generalized circuit theory, we need to make 
a few assumption. First, we neglect spin-flip scattering 
so that the "hat" matrices are purely diagonal (in the 
local basis of the magnetization). Second, we suppose 
that R mx might be non zero but set T mx = 0. Eq. (58) 
and Eq. (59) simplify into 

3La=iRa = ~^T~ 1 ~ VRa] ( 60 ) 

1 _ D L /R 
_1_ " m 
jL/_R,mx - ±^^^/«,- (61) 

where RmJ* are the mixing reflections from left to left 
(R^ x ) and right to right (R* x ). Eq. (60) and Eq. (61) 
are precisely the equations that define generalized circuit 
theory 24 which proves the equivalence with (C)RMT in 
this limit. In fact, the renormalization coefficients of gen- 
eralized circuit theory 24 were chosen such that the cal- 
culation of the conductance with RMT and generalized 
circuit theory fully agree with each other. 



III. A MULTISCALE APPROACH TO SPIN 
TRANSPORT. 

CRMT theory is essentially a spin extension of Ohm 
law and can be considered as an non collinear extension 
of the diffusion equation. Hence, the resistance of a sys- 
tem of length L typically decreases as 1/L and the theory 
cannot describe the exponential suppression of the resis- 
tance found in purely quantum effects such as Anderson 
localization or quantum tunneling. On the other hand, 
these effects (especially the latter) can play an impor- 
tant role in actual devices, such as tunneling magneto- 
resistance (TMR 4,5 ) based spin valves. For those systems 
we need to take a completely different approach from the 
CRMT theory described above and derive a full quantum 
mechanical approach. Such an approach is a priori very 
tempting as it allows to deal with many different systems 
(metals, tunnel junctions, magnetic semi-conductors...) 
within one unified framework and using only few approx- 
imations. However, for practical calculations, such a path 
has a prohibitive numerical cost so that only very small 
systems can be studied (even for a minimum model such 
as the one we will introduce in Section IV). 

Hence, instead of a full quantum general treatment, we 
have a slightly less ambitious double goal: 

(i) we want to study the regime of validity of CRMT 
as well as deviations coming from quantum effects. 

(ii) we want a multiscale description of the system: we 
will provide a purely quantum description of the pieces 
that require it (for instance the insulating oxide layers) 
while the rest of the system is described by CRMT. 
The purely quantum calculation will provide effective 
boundary conditions for CRMT that account for what- 
ever quantum effect is going on there. 

A multiscale approach is rather natural to implement 
in the CRMT framework. Indeed, as the basic ingre- 
dient of the theory are the "hat" matrix, we only need 
to provide the "hat" matrix corresponding to the purely 
quantum parts of the system. Those can in principle be 
calculated through their microscopic definitions, Eq. (6). 
This approach has been applied with good success with 
ab-initio calculations of FN interfaces 30-32 ' 67 . 

In this section, we first show how the "hat" matri- 
ces can be expressed through the intrinsic microscopic 
properties of the system encoded in the retarded Green 
function. Second, we introduce a toy model for a 
ferromagnetic-insulating-ferromagnetic spin valve to pro- 
vide a practical example of the multiscale approach and 
study its regime of applicability. Last, we discuss an ef- 
fective TB model which captures the main features of 
diffusive ferromagnetic layers and compare it with pure 
CRMT calculations. Technically, all the quantum nu- 
merical calculations presented below were done using 
the KNIT package 68 ' 69 within the non equilibrium Green 
function (NEGF) formalism 70 ' 71 . We refer to Ref. 68,69 
for technical details on the numerical algorithm. 
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A. Connection between Semiclassical RMT and 
Quantum NEGF 

Within the NEGF formalism, the transport properties 
of a system are related to the Green function G of the 
system which can be calculated numerically recursively. 
On the other side the "hat" matrices are obtained from 
the microscopic quantum scattering matrix S, as given 
in Eq. (6). In order to connect this two approaches we 
make use of the Fisher-Lee formula 72 which relates the 
scattering matrix S to the retarded Green function G 
and obtain an explicit formula for S in term of the Green 
functions of the system. 

According to Fischer-Lee formula 72 , the transmission 
(reflection) amplitude from a mode m in the the spin 
state a' on the right lead (region 2 in Fig. 1) to a mode n 
in the spin state a on the left (right) lead ( respectively 
region or 2 in Fig. 1) reads, 

= -i^E^G.i^iX^i, (62) 

ij 

ij 

where v n is the velocity of an electron in channels n. The 
summation is taken over interface sites ij which are on 
the border between the system and a given lead, x™j 
(x™j) is the transverse wave function of the m-th mode 
with spin a' evaluated at site j on the Right (Left) lead 
and finally G a i. a 'j is the retarded Green function be- 
tween the spin state a' at site j and spin a at site i. 
Note that the size of the transmission/reflection matrices 
2-^ch x 2iV c h is smaller than the total number of trans- 
verse sites, see Ref. 68 for the definition of the retarded 
Green functions G. 

With these notations, The imaginary part of the reser- 
voirs self energy r L for the Left lead reads, 

K'^kl = J2 hv ^'kX n J- (64) 

n 

with a similar definition for T R . 

The "hat" matrices t is obtained by inserting 
Eqs. (62,64) into Eq. (7), that yields, 

where matrix notation has been used for the summa- 
tion over the site indices. Following the same path with 
Eq. (63), the "hat" matrices f read, 

Tart^'r,' = <W £7777' + ^7~^TrjV ch (T^vGcto-' T^G^, ^ 

+ ^t^c (^r«,G;,,-^r«,G^) (66) 

Using Eqs. (65,66) we can now translate the results 
from a quantum Tight-Binding calculation into our 



CRMT framework and embed quantum TB calculations 
inside CRMT. Alternatively, it will also allow us to com- 
pare pure CRMT with TB calculations and map the 
experimental Valet-Fert parameters Z s f and £j_,Z m x 

into our minimum Tight-Binding model for diffusive 
metallic layers (see Section IV A). 



B. Toy model of a TMR spin valve . 

This section is devoted to studying the accuracy of the 
semi-classical approximation done in RMT and hence of 
the multi-scale approach advocated in this paper. 

The cornerstone of RMT is the "hat" matrix addi- 
tion law for two systems put in serie, as defined by 
Eqs. (16,17). Its validity is based on an ergodic assump- 
tion (an electron entering the system in one channel will 
leave it in arbitrary one) and secondly on the system hav- 
ing many channels so that most quantum effects (such 
as weak localization and universal conductance fluctua- 
tions) can be safely neglected. In term of length scale, 
these two restrictions are related to the mean free path I 
and the localization length £. On one side the size of our 
system L must be long enough (L » £) to reach the dif- 
fusive regime, and therefore suppress the ballistic effects 
that reduce the mixing between channels. On the other 
side the system's size should be small compared to the 
localization length (L <C £) in order to avoid strong local- 
ization, that is not captured by our semiclassical theory. 
Within the TB framework, the inequality £ <C L <C £ 
translates into a window range for the strength W of the 
onsite disorder seen by the electrons. Indeed for the quasi 
one dimensional geometry considered here, l/£ oc W 2 
and £/£ oc N c h. By increasing sufficiently the number of 
open channels, we can be sure to satisfy both require- 
ments. 

If we now add to the system a tunnel barrier, most of 
the open channel will endure a strong reduction of their 
transmission's probability. This corresponds in practice 
to an effective reduction of number of channel really in- 
volved. This render the system more sensible to the influ- 
ence interference effects (universal conductance fluctua- 
tions) and we shall see we will need an even larger number 
of channels to ensure a good agreement between RMT 
and TB calculations. In realistic magnetic nanopillars, 
typically 10 4 — 10 5 channels are involved so that RMT 
is very well justified. However, some care is needed in 
the TB calculations as only a small number of channels 
(a few hundreds) can be considered there. To summa- 
rize, the three main limitation on the applicability of the 
addition law of " hat" matrices are : 

(a) Randomization/Mixing of transverse motion, that 
require a strong enough disorder, (L <C £ ). 

(b) Strong Localization, a purely quantum effect that 
will avoid the use of too strong disorder, (L ^> £). 

(c) Quantum Filtering, that reduces the effective num- 
ber of channel and increases the sensibility to the meso- 
scopic fluctuations. 
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L/2-1 3L/2 




FIG. 4: Sketch of the Ferromagnetic-Insulator-Ferromagnetic 
(FIF) spin valve used to study the validity of the multi-scale 
approach. 



We now consider a tunneling spin valve sketched Fig. 4 
and introduce a TB toy model for its description. The 
ferromagnetic (F) and insulating layers (I) of the FIF- 
system is described with a minimalist TB-model with the 
following Hamiltonian: 



H 



<i,j> 



iv iCi 



(67) 



Here the summation is taken over nearest neighbor sites 
< i,j > and I denote the set of sites included into the 
insulator. The electron destruction operator q is a spinor 
with components: Cj = (cj, c\) where cf annihilates an 
electron with spin a on site i. The discretized system is 
composed of 2L x M x M cells on a cubic lattice, where 
the x axis is the direction of growth of the pillar. The 
onsite energy of TB-model is given by the spin dependent 
scattering potential. 



Vi = Vi 





W l 



(68) 



where the electrons with different spin a feel the same 
profile of the disorder defined by random numbers dis- 
tributed according to a uniform distribution with Vi € 
[—0.5,0.5], but with different strength W . In this toy 
model the tunneling region contain an insulator of length 
of 1 site positioned after L/2 sites. The insulating behav- 
ior is described by the replacement of the onsite energy 
by a spin dependent potential U{ 



Ui = 





u± 



Vi, 



(69) 



In the following, we will compare RMT with pure TB cal- 
culations using two separate calculations: (i) a TB cal- 
culation of the full system of size 2L and (ii) We split the 
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FIG. 5: Comparison between RMT and TB for a single non 
magnetic layer. The error E% between RMT and TB (see 
text) is plotted as a function of the disorder strength W for 
one single sample. The system corresponds to L = 20 and 
M = 5 (iVch = 18, square), M = 6 sites [N ch = 24, diamond), 
M = 8 sites (7V ch = 43, up triangle), M = 20 sites (N ch = 259, 
down triangle). Inset: Error E% as a function of the number 
of open channels iV c h for L = M = 20 with a disorder strength 
W = 1.5. 



system into two parts of size L, calculate the hat matri- 
ces of the two sub part using Eqs. (65,66) and recombine 
them using the RMT "hat" addition law. 

Before dealing with the full spin valve system, we need 
to establish the parameters for which the limitations (a,b) 
are fulfilled. To do so, we perform a spin independent 
survey (W 7 = W) in absence of the insulator (Ui = for 
all sites), and determine the disorder's strength range for 
which RMT agrees with the TB model . In this partic- 
ular case the case the "hat" matrix reduce to scalar and 
addition law from transmission simplifies to, 



RMT — 



-MeftTright 
1 — -Rloft-Rright ' 



(70) 



where Trmt is the transmission probability of the full 
system and Xi e f t ( T^ght) is the transmission probability 
of the left (right) half of the system. We now compare 
the RMT result Trmt with the one directly obtained by 
Ttb calculation obtain for the whole system and com- 
pute the relative error in percent between RMT and TB 
predictions: 



E% = 



100 |Trmt — Ttb\ 
oC^rmt + Ttb) 



(71) 



Note that in order to better observe the discrepancy 
due to fluctuations we perform the calculations on a sin- 
gle sample and do not average over disorder. We checked 
that different samples give consistent results. 

At first the main panel of Fig. 5 shows the relative 
error in function of the disorder strength W. We observe 
the effect of the quantum localization, that corresponds, 





w 

FIG. 6: Comparison between TB (symbols) and RMT (full 
line) of the total transmission T as a function of the disorder 
strength W . M = 10 corresponding to 7V c h = 66 open chan- 
nels and L = 20 (circle), L — 40 (square), L — 60 (diamond). 
Inset: Error i?% in percent between CRMT and TB (see text) 
in function of the ratio L/l. Full symbol correspond to the 
data plotted in the main graph. Empty symbols correspond 
to varying L for fixed W = 2.5 (up triangle) , W = 3.5 (right 
triangle) and W = 4.5 (down triangle). We clearly see that 
we reach the strong localization regime for L/l > 10. 



as predicted, to a clear increase of the error with the 
disorder strength. Such rise of the error can be tempered 
with the increase of number of channels, as indicated by 
the data obtained with N c ^ = 259 (down triangle) . The 
E% as a function of iV c h is investigated in the inset of 
Fig. 5 (-/V c h has been varied at fixed M by changing the 
Fermi energy). Despite of the strong fluctuations, this 
survey clearly shows the strong reduction of the error 
with the increase of N c ^. 

Fig. 6 shows the total transmission trough the sample 
as a function of the disorder strength W for three dif- 
ferent system's size. The data show a good agreement 
between the TB model and CRMT for small disorder. 
Such agreement together with the theoretical prediction 
of an algebraic decay of the transmission as a function 
of mean free path Trmt(^) = (1 — L/l)^ 1 , permit us to 
numerically extract the value of £ = a/W 2 , with a = 40 
and rescales all the length by the mean free path the sys- 
tem's size. The inset of Fig. 6 shows the error in function 
of the ratio 2L/£, for various disorder strength and/or 
system's length. This study indicates that the mixing 
constraint is less important than expected, indeed even 
for a 2L/£ < 0.1 the error is almost null. However, the 
strong localization regime is clearly identified for a ratio 
2L/£ > 10, for which the error exceeds 50%. 

We now turn to the full FIF structure. It is composed 
of an insulator connected to two ferromagnets. The first 
ferromagnets has a size of L/2 — 1 sites and has a fixed 
magnetization, the second one is of size 3L/2 with a free 
magnetization, see Fig. 4. We perform two types of calcu- 
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FIG. 7: Comparison between TB (symbols) and RMT (full 
fine) of the total transmission T of a FIF multilayer as a 
function of angle between layers. We note that L — M = 20 
sites that corresponds to 512 channels, the disorder strength 
is W' =0.5 and = 1.5. The insulator is characterized by 
W — 5 and = 20. Data show a good agreement between 
RMT and TB, we note the error E% is alway bellow one 
precent. 

lations for this sample. First: we calculate the tunneling 
magneto-conductance of the whole multilayer using the 
TB-model, dot point on Fig. 7, as a function of the an- 
gle 9 between magnetization of the two F layers. This 
corresponds to a measure between the probe and 2 in 
Fig. 4. Second: we split the system at the distance L 
(Probe 1 in Fig. 4). We then construct the correspond- 
ing "hat" matrices of each sub-parts and note 5* 01 (S 12 ) 
for area between probes and 1 (1 and 2). We then used 
the " hat" addition law to calculate the " hat" matrices of 
the full system (S 01 " + " S 12 ) and plot conductance, solid 
lines on Fig. 7. Comparison between exact TB-model of 
the FIF layer within RMT show a good agreement, in- 
deed the relative error E% for such unique small sample 
is alway less than one percent, this is quite remarkable 
and validate the multi-scale approach. 

IV. EFFECTIVE TIGHT-BINDING 
HAMILTONIAN APPROACH TO SPIN 
TRANSPORT 

In this section, we consider an effective Tight-Binding 
(TB) Hamiltonian parametrized to reproduce the main 
bulk and interface properties of a few diffusive magnetic 
and non magnetic metals. There are several way one 
can derive an TB Hamiltonian for magnetic multilayers. 
For instance, one can use multi-orbital (s,p,d...) models 
where each site of the TB model corresponds to an atom 
and where all the on site energies and hopping amplitudes 
are parametrized to reproduce ab-initio calculations. 73,74 
When one is interested in transport properties however, 
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one is not interested in a model that reproduces the en- 
tire band structure of a material, but only what happens 
around the Fermi level. Here, we take an even more lim- 
ited approach: we construct our model such that it cor- 
rectly reproduces the 5 characteristic length scales that 
were introduced in the previous section [mean free path 
for majority (If) and minority (ZjJ electrons, spin-flip dif- 
fusion length (l s f), transverse penetration length (l±) and 
Larmor precession length (Zl)]- This quantum model can 
hence be viewed as an effective approach valid for diffu- 
sive metals. It is an extension of the toy model of the 
previous section to properly treat interface properties and 
spin-flip scattering. 



A. A minimum TB model 

We consider a rectangular system of volume L 

X X Ly X 

L z which we discretize on M x x M y x M z cells on a cu- 
bic lattice, each cell having a small volume b x x b y x b z . 
We choose the x axis to be the direction of growth of 
the pillar. Each cell corresponds to one site of the TB 
Hamiltonian which reads, 



tion, with 77?. = ^f- = 1. Different pairs of nearest neigh- 
bors have independent values of 77^ and W so repre- 
sent the strength of spin orbit interaction. The parame- 
ters t?,,t^_, W 1 '' 1 and W so which define the TB-model al- 
low to tune the different characteristic lengths that have 
been introduced before. 

Interfaces. In order to complete the model, the above 
bulk model has to extended to deal with interfaces prop- 
erties. Interfaces are represented as hopping elements 
that link the sites of two different materials. In this view, 
an interface acts as a potential barrier between two ad- 
jacent materials, and does not contain any disorder or 
spin orbit interaction. The hamiltonian that describes 
interfaces is thus simply given by 



Hir. 



£4^ 



(75) 



where the indexes i,j belong to the sites of the adjacent 
materials, and the matrix tp-* is given by 



Ant 



tf 

h 



(76) 



H = -J2 



ck 



i l i3 ^3 



J2 c i ViCi 



<*j> 



£ 

<i,3> 



(72) 



The inter-site hopping matrix tij is material dependent 
and fix the band structure. We use nearest neighbors 
hopping with different values for hopping elements along 
the multilayer growth (x) direction (tn) and within the 
plane (y,z) of the layers (t±). Allowing for different t\\ 
and tj_ is advantageous in the simulations: it allows for 
different discretization steps in the x direction and y, z 
plane. Also, by taking t± significantly smaller than tn, 
one ensures that the number of propagating channels is 
maximum and constant over a large window of energy, 
hence reducing mesoscopic effects. On the other hand, 
the model acquires some intrinsic anisotropy not present 
in the original materials, so that the physical meaning of 
the calculations is doubtful outside the quasi-one dimen- 
sional geometry considered here. These hopping elements 
are also spin resolved, 



H,± = 



t 



t 

1,-L 









(73) 



1,-L 



The second term in the Hamiltonian represents the scat- 
tering potential for electrons with different spins and is 
still described by the Eq. (68). The last term tf° intro- 
duces some spin-orbit interaction responsible for a finite 
spin-flip scattering length. Its spin structure respects 
time-reversal symmetry and is given by, 



t so = W 



(74) 



where the elements of the matrix are random complex 
numbers distributed according to the gaussian distribu- 



Numerical study of a bulk layer within the 
TB-model. 



We are now ready to perform a numerical study of our 
TB model. The quantities of interest are the intrinsic 
total resistance Rj (in unit R s h) for spin u, 



R a i 



2i 



■ Tic 



-1, 



(77) 



the current polarization Pf upon sending polarized elec- 
trons along the a =t, I direction, 



Pf = 



Tfo- — T\ 



Tfa + Tlo- ' 



(78) 



the mixing transmission T mx and the spin-flip probabil- 
ity T-j-j,. Fig. 8 and Fig. 9 show these quantities for two 
sets of parameters that correspond respectively to a nor- 
mal (Au) and a magnetic (Py) material as a function of 
the thickness L of the sample. Together with the quan- 
tum calculations (symbols) we also plot the results of the 
CRMT calculations (lines). The latter is equivalent for 
this collincar case to the Valet Fert Ohmic theory. We 
find that our TB effective model reproduces extremely 
well the semi-classical results: the intrinsic resistance is 
indeed Ohmic (R° = p a L) and the polarization reaches 
its asymptotic value fid after an exponential decay con- 
trolled by Z s f. Exact CRMT expressions for Pf(L) are 
cumbersome but amount with good precision to 



P?(L)^f3 d -(f3 d ±l)e- L ^ 



(79) 



which is well verified by the TB quantum data. For 
magnetic material, T mx is also found to decay exponen- 
tially, as expected. Such an agreement between the TB 
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FIG. 8: Comparison between TB (symbols) and CRMT 
(dashed lines) calculations for transport properties of a nor- 
mal metal as a function of the thickness L of the sample. The 
parameters have been chosen to reproduce the properties of 
gold (Au). Left panel: intrinsic resistance R\(L) = R\{L) 
in unit of the Sharvin resistance. Right upper panel: Pj(L) 
(up triangles) and Pf(L) (down triangles). Right lower panel: 
spin flip transmission probability T-^(L). These results have 
been averaged over 6 samples with different realizations of the 
disorder. 




Material 


p* 




1/irf 


Wt 


Wx 


Wso 


% 


Hi 


t± 


Cu 


5 





0.002 


0.45 


0.45 


0.0055 


1 


i 


0.4 


Au 


20 





0.033 


0.7 


0.7 


0.042 


1 


i 


0.4 


Co 


75 


0.46 


0.017 


4.1 


2.05 


0.022 


3 


i 


0.4 


Py 


291 


0.76 


0.182 


3 


4.6 


0.11 


2 


i 


0.4 



TABLE I: Bulk TB parameters for a few metals. Valet Fert 
resistivity p* is measured in units 10 _9 r2m and spin-flip length 
in nm. The discretization length b z is equal to lnm. These 
parameters have been obtained averaging the calculation of 
the resistance over 6 samples with different realizations of the 
disorder. We have used R B h = 2/Q.m 2 . 




FIG. 10: Resistance of the interface as a function of the hop- 
ping t between two adjacent materials for various mismatches 
of the bulk hopping parameters: Symbols represent the in- 
terface resistance (for a given spin) between two materials 
with i|ii c ft = 3 (red squares), iyieft = 2 (green diamonds) and 



|| left 



1 (black circles). The hopping parameter of the right 



material is kept constant t\ 
of the data of the form Ri n t 
A = 0.09, A = 0.14 and A = 



right 



1. The lines are Pade fits 



= A X (i 2 -t|Heftt||right) /* with 

0.29 respectively. 



FIG. 9: Comparison between TB (symbols) and CRMT 
(dashed lines) calculations for transport properties of a mag- 
netic metal as a function of the thickness L of the sample. The 
parameters have been chosen to reproduce the properties of 
permalloy (Py). Left upper panel: intrinsic resistance Rj(L) 
(circles) and R\{L) (squares) in unit of the Sharvin resistance. 
Left lower: mixing transmission r mx (L). Right upper panel: 
Pj(L) (up triangles) and Pf(L) (down triangles). Right lower 
panel: spin flip transmission probability Tf±(L). These results 
have been averaged over 6 samples with different realizations 
of the disorder. 



C. Interface properties of the TB model 

A similar study can be performed for interfaces proper- 
ties. In this case, the quantity of interest is the interface 
resistance Rf nt (in unit R s h) for spin a, 



per 



(80) 



model and CRMT allows to use the experimentally mea- 
sured CRMT (Valet Fert) parameters to tabulate the TB 
model. The corresponding parameters can be found in 
Table I. We note that the polarization saturates towards 
(3d which defers slightly from the Valet-Fert definition 
of (3 [see Eq. (A14)] due to the depolarizing role of the 
Sharvin resistances. This point will be discussed further 
in Appendix A. 



where T aa represents the transmission probability of the 
interface. The interface resistance has been calculated 
as a function of the hopping at the interface, and 
fitted with simple Pade interpolating formulas (Fig. 10). 
Knowing how the interface resistance varies as a function 
of the hopping allows us to use the Valet Fert experimen- 
tally measured parameters to tabulate the TB model for 
the interfaces, as we have done with the bulk. The TB 
parameters calculated for a few materials are shown in 
Table II. 
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Interface 


r* b (10~ 15 {lm 2 ) 


7 


tf 


h 


Au\Co 


0.5 


0.77 


1.257 


2.2 


Cu\Co 


0.51 


0.77 


1.253 


2.21 


Cu\Py 


0.5 


0.7 


0.99 


2.17 


Au\Py 


0.5 


0.77 


1.257 


2.2 


Au\Cu 


0.5 





1.86 


1.86 



TABLE II: TB parameters for a few Normal-Ferromagnetic 
interfaces. We find parameters of TB model using CRMT 
parameters for the interfaces Eq. (49). The discretization 
length b z is equal to lnm. 



D. Comparison between CRMT and TB for a spin 
valve. 




1 2 1 2 



FIG. 11: Comparison between TB (symbols) and CRMT 
(dashed lines) calculations for the angular dependence of the 
spin torque r exerted on the ferromagnet at the right hand 
side of the stack. Left panel: CusPy2oCusPy2oCu5 multi- 
layer (length in tim). Right panel: CuisgCosCuioPysCiii 
which shows a wavy behavior. 



Now that our model is fully tabulated and that we 
have checked that all the individual pieces (bulk and 
interfaces) agree with CRMT theory, we can go ahead 
and perform quantum calculation for entire spin valves. 
The results are presented in Fig. 11 where we calculate 
the spin torque r (defined as the difference of spin cur- 
rent on the two sides of the free ferromagnet 28,75 ' 76 ) as a 
function of the angle 9 between the magnetization of the 
two magnetic layers. Two different stacks are presented 
(lengths in nm): Cu^Py2QCu^Py2aCu^ (left panel) and 
Cui^gCosCuioPysCui (right panel). The first one is 
symmetric and shows a usual r oc sin 9 torque 75 ' 76 . The 
CRMT and TB calculations are in close agreement. For 
the second stack which is asymmetric, we find a good 
agreement of the torques, except at small angle where 
the TB calculation shows a strong deviation from the 
sin 9 behavior. As a result, the TB torque vanishes for 
a finite value 9 — 9* and the corresponding structure is 
so called "wavy" (see Ref. 77-80 for an extensive dis- 
cussion of waviness where neither the parallel nor the 
antiparallel magnetic configuration is stable). The cor- 
responding CRMT calculation is close, but below, to the 
wavy instability threshold 7 ' 29,81 . This small discrepancy 
indicates a weakness of the effective TB approach as its 
physics is fairly sensitive to the choice of TB parameters. 
To proceed further, we consider in Fig. 12 the symmet- 
ric Cu^Py2oCu^,Py2Q{9)Cu5 stack and compute the spin 
current inside the sample as a function of the position 
x along the stack. Once again, we observe a very good 
agreement between the two approaches. However, this 
new calculation points again to the main weakness of the 
TB approach: the error bars in the TB calculation cor- 
respond to the mesoscopic fluctuations upon averaging 
on different samples (typically 10 in those calculations). 
Real spin valves nanopillars typically have s» 10 4 
channels so that mesoscopic fluctuations ps l/iV c h are 
negligible. Our TB calculations are performed with typi- 
cally 50 — 100 channels, hence do show significant sample 
to sample fluctuations so that quantitative correspon- 
dence with CRMT calculations are only obtained after 
averaging over different disorder configurations. 
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FIG. 12: Comparison between TB (symbols) and C-RMT 
(dashed lines) calculations for the polarization of the spin 
current inside the CusPy2oCu5Py2o(9)Cu5 stack. Magneti- 
zation of the second Py layer is rotated around in XZ plain 
while the first one point in the Z direction. Left, middle and 
right panels stand for 6 = 0,7r/2 and 7r respectively. Upper 
(Lower) panels shows the spacial x dependence of the spin 
current j x = J x /I (j z = Jz/I) along the X (Z) direction. 



To conclude the last two sections, we have found a 
very good agreement between the TB approach and the 
CRMT approach even in regimes where it was not re- 
ally expected: the derivation of CRMT assumes that the 
different channels are ergodically mixed, which is only 
achieved when the typical thickness of the layers is large 
compared to their men free path. In that sense, the 
present TB results can be viewed as a proof of the ro- 
bustness of the CRMT approach as good agreement is 
also obtained in fairly transparent regimes. On the other 
hand, as TB calculations are numerically much more de- 
manding than CRMT, practical calculations are best per- 
formed with the latter approach which do not suffer from 
mesoscopic fluctuations. 
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V. CONCLUSION 

In conclusion, we have given an extensive presenta- 
tion of the semi-classical CRMT theory with an empha- 
sis on the different connections with other existing the- 
ories (Generalized circuit theory and Valet-Fert theory). 
While CRMT is semi-classical, the fact that it has been 
directly derived from the (microscopic) quantum scatter- 
ing approach allows for a simple, direct connection to 
other quantum mechanical calculations. We have used 
this connection to study the regime of validity of the 
semi-classical approximation. In particular we studied its 
accuracy when part of the system (the insulating layer) 
is non-Ohmic and must be described at the quantum me- 
chanical level. The results allows to validate a multi-scale 
approach to magnetic multi-layers where the quantum 
calculations provide boundary conditions for the semi- 
classical simulations. A natural extension of this work is 
to generalize CRMT to three dimensions, allowing to deal 
with system with non trivial magnetic textures such as 
domain wall or vortices 82 Upon embedding the resulting 
theory in micromagnetic simulations, one would reach a 
level of description that should be able to describe quan- 
titatively a large class of spintronic devices. 
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Appendix A: Direct integration of CRMT equations 
Eq. (34) and Eq. (35). 

In this appendix, we provide a direct integration of the 
CRMT equations (34) and (35) for the bulk properties 
of a magnetic layer. An alternative (probably simpler) 
way would be to rely on the mapping with Valet-Fert 
diffusion equation. 

We start with the integration of the "hat" reflection 
matrix defined in Eq. (34) as it is decoupled from the 
"hat" transmission matrix Eq. (35). Let us define the 
stationary solution of ^ — to be f . We seek solutions 
of Eq. (34) in the form: f(L) = f + f\(L) which yields, 



-dL = (roA 



A*)f 1 +f 1 (A r f -A t )+f 1 A'r 1 (Al) 



Since we consider here system with an homogeneous 
magnetization, the inner (mixing) part and of the matrix 
are decoupled from the up and down part so that the 
initial 4x4 matrices can be reduced to 2 x 2. The 2x2 
matrix r (reduced version of r ) satisfies, 



A r - A*r - roA* + r AV = 



(A2) 



where the 2x2 matrices A*/ r are defined in Eq. (43). Af- 
ter some algebra we can transform the preceding equation 
into 



where, 



r = 



1 + r 1 + r 



and Y 





(A3) 



(A4) 



Such kind of quadratic matrix equations can be 
solved 83 and yields 



7-0 



l + r^r-^ryjva 



(A5) 



We proceed with the following change of variable, 



fi(L) = e UBh f(L)e 
A(L) = e ^L A r e n B L 



(A6) 
(A7) 



where Ha = r A r - A* and tt B = A r r - A*. This trans- 
formation simplifies Eq. (Al) into 



dr 



dL 



= -Ml), 



(A8) 



that can be trivially integrated. Inserting the solution 
of Eq.(A8) into Eq.(A6), we obtain the solution for the 
refection "hat" matrix f(L) . The latter can be then 
inserted into Eq. (35) that can be straightforwardly in- 
tegrated to deliver t(L) the transmission "hat" matrix . 
The analytical expressions for f(L) and t(L) are rather 
cumbersome. In what follows we focus on the limiting 
case of L » l s{ where compact expressions can be ob- 
tained. 

Introducing the resistivity per spin channel p^^ as the 
total resistance per spin and per unit length, 



1 



Pa _ j_ 

TZ sh ~ L \T u (L)+T la (L) 



(A9) 



we define the analogs of the Valet Fert parameters for 
the resistivity p* d and spin asymmetry (id as: 



P*d = (Pt + Pi)/ 4 

A, = ^ 



(A10) 
(All) 



Similarly we can determine the relaxation rate of the po- 
larization of the spin current, (it relaxes as e~ L// ' af ) and 
define the corresponding spin- flip length 



Id = -T l0 S 



sf 



T n (L) - T^{L) 



L b \T n (L) + T^(L) 



Pd (A12) 



Re-expressing now Eqs. (A10,A11,A12) with the help 
of our CRMT parameters defined in Eq. (36), we can 
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FIG. 13: Polarization of the spin current inside a Py layer of 
length 150 nm. In the middle of the sample the polarization 
of the spin current is equal to the bulk value /3 (/? = 0.76, red 
dashed line) while close to the ends it decreases to the value 
defined by Eq. (A17), {/3 d = 0.58, black doted line). 



Remarkably only one parameter differ, indeed the spin 
asymmetry /3d contains an additional spin flip length 1 / 
with respect to j3. That difference becomes important 
only when the average resistance within a spin flip length 
is smaller that the Sharvin resistance p*l s i <C TZ s h- This 
discrepancy is due to the implicit presence of the reser- 
voirs in the present calculation which translate (in the 
VF language) into additional Sharvin resistance in series 
with the magnetic layer. This additional Sharvin resis- 
tance induce a depolarization of the current similar to the 
one observed in VF theory at a Ferromagnetic-Normal 
interface where the polarization is reduced to /3' , 



P' = wt, — m ( A16 ) 

with — PnI s { being the resistivity of a normal layer 
within spin flip length. Here we find in analogy, 



directly compare the diffusive values l^,p% Pd to the R 

one obtained with VF-CRMT mapping, as given in 

Eqs. (46,47,48) L + K sh /{*p F l s{ ) 

— = 2y / T\f y/T^ + Tj, = — (A-13) To illustrate this point, Fig. 13 shows the polarization 

^ s f ^ sf of the current inside a Permalloy sample: we find the 

— r-f- . polarization to j5 far away from the boundary and (3d close 

' d " + T_i + 1/1% to t ne reservoirs. This effect should be kept in mind when 

* * comparing quantum results from the Landauer-Buttikcr 

= (r t + rj.)/4. = (A15) approach with Valet- Fert solutions. 
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